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ABSTRACT 

We present new radio continuum data at 4 frequencies on the supermassive, peculiar 
galaxy NGC 1961. These observations allow us to separate the thermal and the non- 
thermal radio emission and to determine the nonthermal spectral index distribution. 
This spectral index distribution in the galactic disk is unusual: at the maxima of the 
radio emission the synchrotron spectrum is very steep, indicating aged cosmic ray 
electrons. Away from the maxima the spectrum is much flatter. The steep spectrum 
of the synchrotron emission at the maxima indicates that a strong decline of the star 
formation rate has taken place at these sites. The extended radio emission is a sign of 
recent cosmic ray acceleration, probably by recent star formation. We suggest that a 
violent event in the past, most likely a merger or a collision with an intergalactic gas 
cloud, has caused the various unusual features of the galaxy. 
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1 INTRODUCTION 

NGC 1961 (Arp 184) is a very massive (dynamical mass 
10 12 - 1O 13 M , Gottesman et al. 1983) Sb spiral galaxy 
with a peculiar and asymmetric optical appearance. It is 
among the most massive and largest known spiral galaxies 
(Romanishin 1983); its optical diameter of D25 = 4.6 ar- 
cmin (de Vaucouleurs et al. 1991, RC3) translates to 110 
kpc (adopting a distance of 82.9 Mpc for Ho = 50 km s~ x 
Mpc -1 ). It is a member of a group of 6 galaxies where it is 
by far the largest and most massive member. 

Despite its size and peculiar appearance, NGC 1961 has 
received relatively little attention. A detailed study of its H I 
structure and kinematics by Shostak et al. (1982, SHSVV) 
found an unusual asymmetric morphology with an extensive 
"wing" of gas extending 30 kpc to the north-west and a sharp 
edge to the south-east. SHSVV suggested that the disturbed 
appearance is due to stripping of the gas in NGC 1961 by 
a hot intergalactic medium (IGM) . They discarded the pos- 
sibility of tidal interaction with another galaxy since there 
is no galaxy with a large enough mass sufficiently close to 
NGC 1961. They also dismissed the possibility of a merger 
with another galaxy or an intergalactic gas cloud since they 
found no trace of a merger remnant such as a second nu- 
cleus, or disturbed velocity distribution in the Hi gas. In 
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favour of the stripping hypothesis SHSVV noted the follow- 
ing points, i) The displacement of the H I with respect to the 
optical disk in the south of the galaxy, ii) The presence of 
increased radio continuum emission in the south-east which 
they interpreted as being caused by an enhancement of the 
magnetic field or by an increase in the cosmic ray (CR) den- 
sity due to a high star formation rate (SFR). iii) The blue 
colour of the optical emission of the NW wing and the south- 
eastern ridge indicating recent star formation, iv) A curved 
X-ray emission feature in the south-east of the galaxy de- 
tected by the Einstein satellite (Shostak et al. 1982) which 
they interpreted as emission from hot IGM. However, recent 
ROSAT observations have shown that the X-ray emission of 
NGC 1961 is fairly typical for a spiral galaxy and comes 
mainly from the central region (Pence & Rots 1997). The 
cause of the asymmetric Hi distribution therefore remains 
unclear. 

In this paper we present multi-frequency radio contin- 
uum data which we use to identify thermal and nonther- 
mal (synchrotron) emission. The thermal radio emission re- 
veals the sites and intensity of present day star formation, 
whereas the synchroton emission, and in particular its spec- 
tral shape, contains information about the star formation in 
the past (on time-scales of a few 10 7 yr which is related to 
the life-time of the CR electrons and of supernova (SN) pro- 
genitors), and also about the acceleration and propagation 
mechanisms of CRs. 
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2 OBSERVATIONS 

NGC 1961 was mapped with the VLA| at 8.41 GHz and 
with the Cambridge 5-km Ryle Telescope (RT) at 4.92 and 
15.4 GHz. Additionally, existing VLA data (B-array) at 
1.46 GHz (Condon 1983) were available. When the obser- 
vations ol NGC 1961 were made (between 1990 and 1993), 
the RT had a band-width ol 280 MHz split into 28 10-MHz 
frequency channels; together with the minimum baseline of 
18 m this provides excellent temperature sensitivity at both 
4.92 and 15.4 GHz. 

For the RT observations a total of five telescope con- 
figurations at each observing frequency were employed giv- 
ing a nearly fully filled aperture out to a baseline of ap- 
proximately 1 km. A minimum of two 12-hour runs in 
each configuration were obtained. This gave a resolution of 
approximately 1.0 x 1.0cosec(<5) arcsec 2 at 15.4 GHz and 
3.0 x 3.0cosec(<5) arcsec 2 at 4.92 GHz. At 4.92 GHz observa- 
tions of phase-calibrators were made at the beginning and 
end of each run, while at 15.4 GHz calibration observations 
were interleaved with those of NGC 1961. 3C 286 and 3C 48 
were observed regularly as flux calibrators. Calibration and 
data-editing was performed in the MRAO package post- 
mortem, with subsequent reduction in MPS and the MRAO 
package ANMAP. 

The 8.41-GHz VLA observations were in the B, C and 
D arrays and were carried out during 1990/91. Reduction 
of the data followed standard VLA procedures with calibra- 
tion, editing and imaging performed in AIPS. The 1.46 GHz 
data were from Condon (1983) and no additional reduction 
was required. A number of maps were made at each fre- 
quency by tapering the aperture-plane to highlight structure 
on different angular scales. 

The data in this paper which are most seriously affected 
by the lack of short baselines are those at 1.46 GHz, which 
are only sensitive to structures on scales smaller than about 
2 arcmin. Very faint emission from the whole of the optical 
disc of NGC 1961 will not be adequately represented in our 
maps. The optical diameter of NGC 1961 is 4.6 arcmin; how- 
ever the extent of the continuum emission at 1.46 GHz as 
determined by Shostak et al. (1982) using the WSRT with a 
shortest baseline of only 36 m is approximately 2.5 arcmin, 
with very little low-brightness emission outside that region. 
A further test of the completeness of our mapping is to com- 
pare our measured flux densities with single dish data; these 
comparisons are reported in Table 1. It is apparent that to 
within the quoted uncertainties our flux densities are con- 
sistent with those determined by other workers (including 
single dish measurements at 1.40 GHz and 4.85 GHz), al- 
though the best estimates represent a shortfall of 18% of 
the flux density at 1.46 GHz and 12% at 4.92 GHz. These 
results suggest that we are able to image the extent of the 
continuum emission satisfactorily at all our observing fre- 
quencies. 

So as to make the spectral index comparison as accu- 
rate as possible it is crucial to ensure that the spectral index 
is determined from images all of which are sensitive as far as 
possible to the same range of spatial structure. To achieve 

t The VLA is operated by the National Radio Astronomy Ob- 
servatory for Associated Universities Inc., under a cooperative 
agreement with the National Science Foundation. 



Table 1. Integrated flux densities 



Frequency Flux Density Reference 



IIVLXIZI 






38 


3.40 ± 0.85 


Howarth (1990) 


57.5 


2.0 ± 0.5 


Israel & Mahoney (1990) 


151 


1.17 ± 0.21 


Howarth (1990) 


1400 


0.191 


White & Becker (1992) 


1465 


0.149 ±0.015 


this work 


1490 


0.182 ± 0.023 


Condon (1987) 


4850 


0.057 ± 0.009 


Becker, White & Edwards (1991) 


4919 


0.050 ± 0.006 


this work 


8414 


0.028 ± 0.003 


this work 


10550 


0.031 ± 0.005 


Niklas et al. (1995) 


15360 


0.019 ± 0.002 


this work 



such a match, the UV-range of the data in the aperture plane 
used in the construction of the images above 1.46 GHz was 
restricted during the mapping programs so that the min- 
imum and maximum baselines (measured in wavelengths) 
used in the imaging were the same; this involved removing 
baselines from all of the high-frequency data. While it is not 
possible to ensure that the precise sampling of the aperture 
plane is the same in all cases, for both the VLA and RT the 
dense sampling near the shortest baselines leads to similar 
overall coverage of the UV plane. Furthermore the use of 
CLEAN does, to some extent, compensate for the resulting 
marginally different beamshapes. A possible problem results 
from the fact that the extended, low-brightness structure of- 
ten has a steep radio-frequency spectrum. Missing the short 
baselines could lead to a negative region ("negative bowl") 
around the source and hence a (frequency-dependent) error 
in the zero level of the image. The final images were checked 
for the existence of such an effect and in each case any offset 
of the local zero level was found to be less than the thermal 
noise (this is clear from Figs. 1 and 2 where no evidence of a 
significant negative bowl exists). We are confident that the 
procedures that we have adopted here are sufficient to reduce 
any artefacts resulting from the sampling of the aperture to 
levels significantly below those of the thermal noise. 

All the maps were CLEANed and a set of maps at 
a common resolution of 16 x 16 arcsec 2 were used for the 
mult i- frequency comparison (at 1.46 GHz this was achieved 
by convolving the supplied image) . All maps have been cor- 
rected for the primary beam response of the telescope used 
to make the observations. 



3 RESULTS 

3.1 Radio images 

Fig. 1 shows the radio maps at the four frequencies at a 
resolution of 16 arcsec; maps at a resolution of 8 arcsec are 
shown in Fig. 2. In Fig. 3 the 1.46 GHz radio map is overlaid 
on an optical image obtained from the Digitized Sky Survey 
produced by the Space Telescope Science Institute. The ra- 
dio images at all frequencies show a peculiar morphology 
which are similar to the optical image. The bright central 
nucleus is resolved in the higher resolution images and con- 
sists of two peaks - one which corresponds to the optical 
peak and a second peak which coincides with a dust lane. 
Furthermore, a bright radio arm is visible in the south-east, 
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Figure 1. Maps of NGC 1961 at 1.46, 4.92, 8.41 and 15.4 GHz at a resolution of 16 X 16 arcscc 2 . The contour levels start at (-0.85, 

-0.60, 0.60, 0.85 ) mjy beam- 1 (1.46 GHz and 4.92 GHz); (-0.28, -0.20, 0.20, 0.28,...) mjy beam- 1 (8.41 GHz) and (-0.48, -0.34, 0.34, 

0.48....) mjy beam" 1 (15.4 GHz) and increase by factors of \/2. The rms noise levels are: 0.30 mjy beam -1 (1.46 and 4.92 GHz), 0.10 
mjy beam- 1 (8.41 GHz) and 0.17 mjy beam- 1 (15.4 GHz). 

1.46 GHz 8.41 GHz 
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Figure 2. Maps of NGC 1961 at 1.46 and 8.41 at a resolution of 8 X 8 arcscc 2 . The contour levels start at (-0.48, -0.34, 0.34, 0.48 ....) 

mjy beam" 1 (1.46 GHz) and at (-0.14, -0.10, 0.10, 0.14 ) mjy beam" 1 (8.41 GHz) and increase by factors of \/2. The rms noise levels 

arc: 0.17 mjy beam- 1 (1.46 GHz), 0.05 mjy beam- 1 (8.41 GHz). 
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Figure 3. The radio map at 1.46 GHz from Fig. 2 overlaid on an 
optical image. The contour lines start at 0.34 mjy and increase 
by factors of 2. 

coinciding with the bright optical south-eastern arm. The 
optical spiral arm extending from of the west of the nucleus 
and to the south is also visible in the radio image. At all 
frequencies, except for 15.4 GHz, there is an extended enve- 
lope of radio emission visible which surrounds the nucleus 
and the south-eastern arm and which contains, at 1.46 GHz, 
about 50% of the total radio flux. 

3.2 Radio spectrum and spectral index 
distribution 

Integrated flux densities at the four observing frequencies 
were determined from the low-resolution maps and in Table 
1 the results are summarized together with flux densities 
taken from the literature. The integrated spectrum is shown 
in Fig. 4; the spectrum is well-fitted by a power law with a 
spectral index a — 0.85 (S(u) oc v~ a ), which is not unusual 
for a spiral galaxy; the average spectral index for a sample 
of spiral galaxies is 0.74 ±0.03 (Gioia et al. 1982). 

As discussed in Section 2 there may be some flux miss- 
ing from our images due to the restricted UV-coverage on 
the shortest baselines; our estimate of this is 18% (33 mjy) 
at 1.46 GHz and 12% (7 mjy) at 4.92 GHz. The spectral 
comparison has been performed using images with matched 
UV-coverage and resolution so that the spectra of structure 
of similar spatial scales may be compared. If the emission is 
regarded as features superimposed on a background of emis- 
sion, then on those scales which are properly represented in 
our imaging the spectral results will be accurate. A prob- 
lem may exist, however, in those regions which do not show 
significant structural variation. In this case we can quantify 
the likely errors. If we take a worst case of 33 mjy miss- 
ing flux density at 1.46 GHz distributed over an area of 9 
armin 2 (slightly larger than the observed mappable extent 
of the source), then this would correspond to an error of 
0.26 mjy/beam over the source at 1.46 GHz, which is less 
than the thermal noise and at a level of approximately 40% 
of the lowest contour in Fig. 1. Since we would expect the 
spectrum of the most extended material to be steep, this 
effect becomes less important at the higher frequencies. The 




log(i/) [GHz] 

Figure 4. The radio spectrum of NGC 1961. The data are listed 
in Table 1. The line is the best- fit power-law spectrum and has a 
slope of 0.85. 



effect of the missing short spacings is therefore in the worst 
case to underestimate very slightly the flux density at 1.46 
GHz relative to the higher frequencies and therefore to un- 
derestimate also the spectral index. However, this effect will 
lead to a systematic error which is considerably smaller than 
the statistical error due to the thermal noise and in what fol- 
lows it will be ignored. 

Fig. 5 shows spectral index maps between 1.46 and 4.92 
GHz and between 1.46 and 8.41 GHz, each at 16-arcsec reso- 
lution and gated so that the spectral index is only calculated 
where the flux density exceeds 3<r at both frequencies. The 
spectral index between 1.46 and 4.92 GHz (ai.46-4.92GHz) 
shows a peculiar behaviour in the sense that at the max- 
ima of the radio emission (i.e. at the nucleus and on the 
south-eastern arm) the spectral index is steep. Away from 
these sites the spectral index flattens, except towards the 
north and west of the nucleus and towards the north of 
the arm. The behaviour of ai.46-8.4iGHz is less uniform 
which is largely due to the contribution of thermal emis- 
sion at the higher frequency. The most noticeable flattening 
of the spectral index in Fig. 5 is between the nucleus and 
the south-eastern arm where Qi.46-4.92GHz = 0.2 ± 0.3 and 
Qi.46-8.4iGHz = 0.55±0.2. The spectral index variation that 
would normally be expected due to synchrotron and inverse 
Compton aging of the CR electrons is a gradual steepening 
of the spectrum away from the sites of acceleration. How- 
ever, at the distance of NGC 1961 our 16-arcsec beam cor- 
responds to 6.4 kpc which is larger than the typical scale 
of electron diffusion, and we are effectively averaging over a 
CR electron population of varying ages. 
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Figure 5. Maps of the spectral index between 1.46 and 4.92 GHz (left) and between 1.46 and 8.41 GHz (right) at a resolution of 16 X 16 
arcsec 2 . The greyscales range from 0.1 to 1.2 (ai. 46-4. 92GHz) an d from 0.6 to 1.2 (ai.46-8.4lG.Hz)- Overlaid is a 1.46 GHz image with 
contour lines starting at 0.6 mjy and increasing by factors of 2. 



4 SEPARATING THE THERMAL AND 
SYNCHROTRON EMISSION 

The maps at 4 frequencies allow us to separate the radio- 
continuum emission into thermal and nonthermal compo- 
nents. A spectral-fitting technique is used to estimate the 
two contributions and determine spatial variations in the 
form of the synchrotron spectrum. For a set of maps at a 
common resolution a function of the following form is fitted 
at each pixel 



I(u) = SC(x B ,j)+T(iy )iy- 



(1) 



where T(z^o)^ 01 is the emission from optically thin ther- 
mal Bremsstrahlung and SC(xb,j) is the synchrotron emis- 
sion from an aged electron population with injection index 
7 which has been radiating for a time t; xb = p/vb, where 
vb is the break frequency. During the fitting procedure the 
minimum in \ 2 is found by varying S, vb and T(vq) with 
the additional constraint that they are positive definite. The 
form of the synchrotron spectrum C(xb, 7) depends on many 
processes, but to simplify the problem we make the following 
approximations. 

(i) The dominant electron energy loss mechanisms are 
synchrotron and inverse Compton emission. These energy 
losses are proportional to the energy density of the mag- 
netic field, Ub, and the energy density of the radiation 
field, C/ ra d, respectively. The magnetic field can be estimated 
from the standard minimum energy assumption. With a 
high-frequency cut-off of 10 GHz, a low-frequency cut-off of 
100 MHz, a spectral index of 0.8, and a ratio of CR proton- 
to-electron energy density of k — 100, this yields B = 15/iG. 
[Trad can be estimated from the bolometric luminosity, Lboi, 
of the galaxy by assuming that the radiation emitting mat- 
ter is distributed in a flat cylinder: 



2L h 



1R 2 



(WW) 35 

— L x 1.4 x 10" 35 

(R/kpc) 



eV 



(2) 



We approximate L boi as the sum of the far-infrared and 
blue emission (yielding L b oi = 1.3 1 38 W) and hence find 
f/rad = 1.5 eV cm~ 3 . 



(ii) Diffusion in the disc of the galaxy is negligible on 
scales larger than our observing beam, which in this case is 
a good approximation. 

(iii) The electrons which are radiating within a given res- 
olution element were all accelerated (injected) a time t ago 
(the electron age). 

(iv) The electron loss time is much greater than the pitch 
angle randomization time — the so-called JP model (Jaffe 
& Perola 1974). 

(v) The electron injection index is 7 = 2.0. 

With these assumptions the electron energy distribution 
with an injection index 7 is given by: 



N(E,t) = 



N E-~<(1 - e(U B + U^Ety- 2 E<E B 
E>E B 



(3) 



with e = (4/3)(<JTc/mc 2 ), ar being the Thompson cross 
section, and Eb = {(.(Ub + t/ ra d)t) _1 . The synchrotron spec- 
trum is obtained by convolving this electron energy distri- 
bution with the spectrum for a single electron. The break 
frequency is given by: 



vb = (3/2)v G E B oc 



(f/rad + U B )H 2 



(4) 



where vq is the gyrofrequency. The fitting procedure was 
applied to all pixels at which the flux density exceeded 3 
times the noise at 1.46, 4.92 and 8.41 GHz and 1.5 the noise 
at 15.4 GHz. The lower gate value at 15.4 GHz was chosen 
in order included the extended, predominantly nonthermal, 
radio emission in the analysis. Since the fitting procedure 
takes full account of the error in the data, the quality of 
the separation is not affected by this. In Fig. 6 the results 
of this separation are shown in the form of maps of the 
thermal emission, synchrotron emission and electron age de- 
rived from the breaks in the synchrotron emission. The syn- 
chrotron emission is more extended than the thermal emis- 
sion which is mainly concentrated around the maximum of 
the radio emission and the south-eastern arm. The electron 
age reflects the peculiar behaviour which is seen in the spec- 
tral index maps and noted above. The ages are highest at 
the maxima of the radio emission where the spectrum be- 
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Figure 6. Maps of the synchrotron emission (left), thermal emission (middle) and break frequency (right) from which the age of the 
electrons can be derived. For the synchrotron emission the greyscale ranges from (white) to 19 mjy/beam (black), for the thermal 
emission from (white) to 4 mjy/beam (black) and for the break frequency from 1 (white) to 20 GHz (black) (corresponding to electron 
ages from 1.4 10 7 yr (white) to 3 10 6 yr (black)). 



tween 1.46 and 4.92 GHz is steepest. The total thermal radio 
continuum emission found from this fitting procedure is 16 
mjy at 1 GHz (i.e. ~ 10% of the total emission). 

Fig. 7 shows selected radio-continuum spectra of 
NGC 1961 and the fits associated to these data together 
with their positions relative to a 1.46 GHz map. The steep 
spectrum at the nucleus and at the maximum of the south- 
eastern spiral arm can be seen (positions 3 and 4) . At these 
locations the flattening of the spectra towards higher fre- 
quencies due to a high thermal component is clearly visi- 
ble and the thermal emission can be determined with good 
accuracy. The thermal fraction is relatively high (25% at 
1 GHz) and the nonthermal spectrum steep (a n th = 1-7 be- 
tween 1.46 and 4.92 at both positions). The derived break 
frequency is low, suggesting that acceleration has not taken 
place for a considerable period. Away from these positions, 
the radio spectra show less flattening towards higher fre- 
quencies (positions 5 and 6), or none at all (position 1 and 
2) and therefore indicate a lower thermal contribution. The 
radio spectra at the positions 1 and 2 are flat at low fre- 
quencies and steepen toward higher frequencies. The break 
frequency is high, indicating recent acceleration of CR elec- 
trons. 



5 DISCUSSION 

5.1 FIR and radio luminosity 

In Table 2 observational properties of NGC 1961 are pre- 
sented. Because of its size, NGC 1961 is luminous at all 
wavelengths. The {B — V)% colour is typical for an Sb galaxy 
(Huchra 1977) as is the ratio Lfir/£b = 0.3 (e.g. Lisenfeld 
ct al. 1996b). This latter value is a measure of the ratio of 
the SFR averaged over the last 10 8 yr to the SFR averaged 
over the last 10 9 yr, and indicates that the global SF has 
proceeded constantly during the past 10 9 yr. 

The ratio of the far-infrared (FIR) to radio luminosity, 
log(LFm/ PiA9GHz) = 14.44, is a factor of 3-5 lower than 
for spiral galaxies. The average values derived for galaxy 
samples range from 14.9 (bright galaxies, Condon Ander- 
son & Helou 1991) to 15.1 (normal spiral galaxies, e.g. 
Condon et al. 1991, Lisenfeld et al. 1996a) with a typi- 
cal standard deviation of 0.2 in the logarithm. This result 
is most likely due to an unusually high radio luminosity; 
the FIR luminosity does not appear unusually low com- 
pared to its blue luminosity (Lfir/£b = 0.3), or the in- 



Table 2. Global properties 



Distance 


82.9 Mpc 


SHSVV 1 


£>25 


4.6 arcmin 


RC3 


Inclination 


50° 


SHSVV 


-Pl.49G.Hz 


1.5 10 23 W Hz" 1 


Condon (1987) 


SlOO 


6.4 Jy 


IRAS Faint Source Cat. 


SeO 


21.6 Jy 


IRAS Faint Source Cat. 


Seo/Sioo 


0.30 




Lfib. 2 


3.9 10 37 W 




B°, 


11.01 


RC3 


Lb 


9.5 10 37 W 




(B - V)° 


0.57 


RC3 


M Hl 


1.2 10 11 


SHSVV 1 


M Ha 


4 10 10 


Young et al. (1995) 3 


M tota l 


2.3 10 12 


SHSVV 1 



1 adjusted to Hq = 50 km s 1 Mpc 1 

2 L FIR = 1.5 10 32 [D/Mpc) 2 (2.58(S 60 /Jy) + (Sioo/Jy))] 

3 assuming a conversion factor of Mjj n /Leo = 4.8 

ferred H2 mass (L fir /Mh 2 = 2LQM0 1 ; for our Galaxy 
L Fm /M H2 = 3.7LQMQ 1 , Solomon et al. 1992). On the other 
hand, the radio luminosity is high compared to luminosities 
in other bands. For example the radio to blue luminosity ra- 
tio P rad /L B = 6.25 lO'^Hz-'L' 1 is a factor of 7.5 higher 
than found for a sample of Sb galaxies (Hummel et al. 1988). 
Furthermore, the radio surface brightness is high as reflected 
in the magnetic field strength, B = 15piG, estimated using 
the minimum energy method in which B oc (P ra( j/area) 2 ^ 7 . 
Hummel et al. (1988) derived for a sample of Sb galaxies 
(with the same assumption for the minimum energy calcu- 
lation) an average value of 8 fiG. 

Deviations from the FIR/radio correlation have been 
found before for cluster (Niklas, Klein & Wielebinski 1995) 
and interacting galaxies (NGC 2276, Hummel & Beck 1995). 
Also in these cases the conclusion by the authors was that 
the radio emission is enhanced compared to normal galax- 
ies. The reasons for this enhancement have been suggested 
to be either an additional production of CRs by shock accel- 
eration at interstellar shocks produced by the interaction of 
the galaxy with the intercluster medium (Volk & Xu 1993), 
or an enhanced magnetic field (Hummel & Beck 1995). 



5.2 The spectral analysis 

The diffusion scale of the CR electrons in NGC 1961 is 
approximately 3.7 kpc at 1.46 GHz (taking B — 15^iG, 
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Figure 7. Selected spectra in NGC 1961 together with the positions where they were extracted relative to a 1.46 GHz image. 



t/rad = 1.5 eV cm , and a typical diffusion coefficient of 
10 29 cm 2 s _1 ). At the distance of NGC 1961 the 16-arcsec 
beam corresponds to 6.4 kpc; to a good approximation, 
therefore, the observations average over the synchrotron 
emission of both newly accelerated and aged CR electrons. 
The spectral variations reflect variations in the thermal frac- 
tion, differences in the conditions of the ISM and in the CR 
injection, but not effects due to propagation of the CR elec- 
trons. 



5. 2. 1 Star formation history 

The shape of the non-thermal spectrum reflects temporal 
variations of past CR acceleration and therefore of the star- 
formation history if SN remnants are the most important 
source for the acceleration of CR electrons. The steep syn- 
chrotron spectrum at the centre and in the south-eastern 
arm (positions 3 and 4 in Fig. 7) indicate aged CR elec- 
trons; the spectra cannot be explained by continous CR 
acceleration since this would require an injection spectral 
index of 7 = 3.4. In the region between the centre and the 
south-eastern arm the synchrotron spectrum is much flat- 



ter, indicating recent electron acceleration. The electron age 
distribution derived from the spectral fitting analysis gives 
a time since acceleration for the CR electron population at 
the maximum of the radio emission and in the south-eastern 
arm (positions 3 and 4) of ~ 10 7 yr, whereas the electrons 
between these sites are considerably younger (3 x 10 6 yr). 

This leads to the following picture of the star forma- 
tion history. The end of CR injection 10 7 yr ago indicates 
that the SN activity ceased at that time. The life-time of SN 
progenitor stars is at maximum 8 x 10 7 years (correspond- 
ing to a stellar mass of 5 Mq), which suggests an epoch of 
intense SF which ended <, 9 x 10 7 yr ago. The beginning 
of the star formation epoch cannot be constrained by our 
data. The recent acceleration of CR electrons deduced from 
the flat spectral index of the extended radio emission means 
that the CR injection and thus SN activity at these sites has 
recently started. The spectra are too flat to be explained by 
continous CR acceleration because in this case the contri- 
bution from aged CR electrons to the radio emission would 
steepen the spectrum, and approximating this situation to 
that of a steady-state we would expect a synchrotron spec- 
tral index of a nt h = 7/2. Taking 7 = 2, the expected syn- 
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chrotron spectral index on the basis of a steady state model 
is much steeper than the observed one, which is, attributing 
as an upper limit all of the observed flux density at 15 GHz 
to the thermal emission, a nt h, 1.46-4. 92GHz = 0.4±0.2. These 
results are most consistent with a recent burst of star for- 
mation 2 — 3 x 10 7 yr ago followed by a peak SN rate about 
5 x 10 6 yr ago. 

The time-scales derived from this analysis imply that 
a causal connection between the star formation epochs is 
possible, in the sense that the earlier star formation period 
in the arm may have triggered an episode of star formation 
in the inter-arm region. The width of the south-eastern arm 
is about 10 kpc and taking a typical sound speed to be 10 
km/s a SN shock would take approximately 10 yr to cross 
the arm which is consistent with the time-scale inferred for 
the star formation history. 

SHSVV suggested that an encounter of NGC 1961 with 
an intergalactic cloud took place less than 5 10 8 yr ago (the 
rotation time-scale), and triggered an episode of intense star 
formation in the south-eastern arm. They presented an op- 
tical colour image showing that the south-eastern arm has 
indeed a blue colour, indicating recent star-formation. This 
is consistent with the star-formation history that we infer. 



Aa = a - = 
2 



I f (2 + 7) 
2 2+f 



(6) 



For an increasing magnetic field with a = — 1 this yields 
Aa = 0.4 and for a — —2 Aa = 0.66. Thus a signifi- 
cant steepening of the spectrum is possible. If however syn- 
chrotron losses (oc B 2 ) dominate then, since U ra d + Ub ~ 
Ub oc B 2 , we set p = 2a and get: 



Aa = 



*(1- |(2 + 7)) 



(7) 



In this case only a very weak steepening of the spectrum is 
possible; for 7 = 2 Aa is zero for all values of p and a. 

For NGC 1961 the high radio surface brightness sug- 
gests that synchrotron losses dominate. In this case, the 
above analysis shows that under steady-state conditions 
variations in the field strength do not lead to significant vari- 
ation in the spectral index. Steep spectra in the south-east 
arm could, however, result from a non-steady-state scenario 
in which enhancement of the field was followed by a burst 
of star formation with increased synchrotron loss rate - it 
is not possible to modify the spectrum significantly without 
incorporating into a model variation in the star-formation 
rate. 



5.2.2 Variations of the magnetic field 

In the above analysis we have assumed that the magnetic 
field strength is constant over the galactic disk; this is 
very likely not the case. Compression of the ISM and/or 
the presence of shocks throughout the disc will lead to en- 
hanced magnetic field strengths and hence increased radio 
synchrotron luminosity relative to FIR. Such processes were 
inferred by SHSVV from the morphology of the H I distribu- 
tion. However, an enhanced magnetic field strength in the 
south-eastern arm alone cannot explain the large radio-to- 
FIR ratio since the radio flux density at 1.46 GHz of the 
south-eastern arm contributes only about 15 % to the total 
flux density of the galaxy. 

If the magnetic field varies on scales shorter than that 
of the CR electron propagation it has an effect not only on 
the local synchrotron emission but also on the integrated 
spectral index. In particular, if the magnetic field strength 
decreases (increases) away from the source of CR electrons, 
the spectrum can become significantly flatter (steeper). The 
problem can be solved in a one-dimensional model in which 
the magnetic field strength and the energy density of the ra- 
diation field change with increasing height over the galactic 
disk, z. Assuming that the magnetic field strength varies as 
B(z) oc z~ a , the energy losses as ((7 rat j + Ub){z) oc z~ p 
and the energy dependence of the diffusion coefficient is 
D(E) oc , then the resulting spectral index is, in the 
steady state, 



7 1 - 



H P 



r(2 + 7) 



2-p+f(l- M ) 



(5) 



(Berezinsky et al. 1990). Taking p — 0, we see that in the 
case p — (i.e. spatially constant inverse Compton losses 
dominate the energy losses) the difference of the spectral 
index with respect to the case of a constant magnetic field 
is 



5.3 A dynamical history for NGC 1961 

SHSVV suggest that a collision with an intergalactic cloud 
has resulted in the unusual Hi distribution in NGC 1961. 
Such an interaction with the IGM can explain in addition 
to the optical and H I morphology many features of the radio 
emission provided the interaction triggered an episode of in- 
tense SF in the south-eastern arm. The main difficulty with 
this model is that the massive intergalactic cloud required 
has not been observed either as ionized (Pence & Rots 1997) 
or in neutral (SHSVV) gas. 

An alternative view is that the disturbed appearance of 
NGC 1961 could be the result of a merger. SHSVV discarded 
this possibility since they found no evidence for it either in 
the optical (i.e. there is no second nucleus) or in the H I 
distribution. However, our high resolution radio images re- 
veal the presence of a second nucleus in the centre of the 
galaxy consistent with an advanced merger; the disturbed 
optical and Hi appearance is then naturally explained as a 
merger remnant. Near-infrared observations will provide an 
important indication as to whether a highly obscured optical 
nucleus is present. 



6 CONCLUSIONS 

We have presented radio data at 4 frequencies of the su- 
permassive spiral galaxy NGC 1961. Using these data we 
fitted a combined aged-synchrotron plus thermal emission 
spectrum to each point in the galaxy. The most recent star 
formation, that is traced by the thermal radio emission, is 
mainly taking place in the centre and the south-eastern arm 
of the galaxy, whereas the underlying extended radio emis- 
sion is predominantly nonthermal. 

The spectral index distribution between 1.46 and 
4.92 GHz over the disk of the galaxy shows a peculiar be- 
haviour. At the maxima of the radio emission the nonther- 
mal spectral index is very steep whereas the extended radio 
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emission has a much natter spectrum. We discuss various Volk H.J., Xu C, 1993, Infrared Phys.Technol. 35, 527 

possibilities to explain this peculiar behaviour. The most White R.L., Becker R.H., 1992, ApJS, 79, 331 

likely explanation is that the steep spectrum is due to the Young J., Xie S., Tacconi L., et al., 1995, ApJS, 98, 219 

end of the CR injection about 10 7 years ago. Assuming that 
CRs are accelerated mainly in SNRs this means that a phase 
of intense SF has ceased <, 10 8 yr ago. The flat spectrum 
of the extended emission is a sign of recent CR acceleration 
in these regions. This could be due to CR acceleration in 
the shocks of SNRs of a recently formed stellar population. 
The radio-to-FIR ratio exceeds the average value of spiral 
galaxies by a factor of 3-4. An overall increase of the mag- 
netic field with respect to the radiation field could account 
for this high FIR/radio ratio. Spatial variations of the mag- 
netic field, which are likely to be present if the magnetic 
field is frozen into the gas, are however unable to give the 
observed distribution of the spectral index. 

The discovery of a second radio nucleus leads us to sug- 
gest that NGC 1961 is most likely the result of a merger 
rather than collision with a massive intergalactic gas cloud. 
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